Electronic instabilities in 3D arrays of small-diameter (3,3) carbon nanotubes 
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We investigate the electronic instabilities of the small-diameter (3,3) carbon nanotubes by study- 
ing the low-energy perturbations of the normal Luttinger liquid regime. The bosonization approach 
is adopted to deal exactly with the interactions in the forward-scattering channels, while renormal- 
ization group methods are used to analyze the low-energy instabilities. In this respect, we take 
into account the competition between the effective e-e interaction mediated by phonons and the 
Coulomb interaction in backscattering and Umklapp channels. Moreover, we apply our analysis to 
relevant experimental conditions where the nanotubes are assembled into large three-dimensional 
arrays, which leads to an efficient screening of the Coulomb potential at small momentum-transfer. 
We find that the destabilization of the normal metallic behavior takes place through the onset of 
critical behavior in some of the two charge stiffnesses that characterize the Luttinger liquid state. 
From a physical point of view, this results in either a divergent compressibility or a vanishing renor- 
malized velocity for current excitations at the point of the transition. We observe anyhow that this 
kind of critical behavior occurs without the development of any appreciable sign of superconducting 
correlations. 



perature. 

Nevertheless, superconductivity at about 15 K has 
been claimed to occur in 4 A-diameter nanotubes^. 
In the experiment reported in Ref. 0, a strong dia- 
magnetic behavior was interpreted as an anisotropic 
Meissner effect, while a genuine superconducting tran- 
sition was not observed. This has opened some contro- 
versy on this issue, since ab initio simulations predict a 
room-temperature Peierls transition in the allowed 4A- 
diameter geometries, namely in the (5,0)A^ and the (3,3) 
nanotubes^. On the other hand, mean-field calculations 
for the (5,0) nanotubes seem to find a superconducting 
instability, but with a critical temperature of about 1 
Kit 



The development of nanoscale technology during the 
last decade has attracted much attention on carbon nan- 
otubes, which are among the most promising candidates 
to fabricate molecular-size devices. This is mainly due to 
the wide variety of their electronic and transport prop- 
erties, which can result in metallioi, semiconducting^ or 
even superconducting behavior—, depending on geometry 
and the way of assembling. 

^From a theoretical point of view, the confinement of 
electrons in the longitudinal dimension of the nanotubes 
induces the so-called Luttinger liquid behavioipiiL&Siill. 
This is characterized, for instance, by the power-law de- 
pendence of the differential conductance, which has been 
actually observed experimentally!^. 

Such a behavior breaks down anyhow at sufficiently 
low temperature, and the nanotubes enter a different 
regime, usually driven by the quality of the contacts 
in the experimental setup. In particular, in the case 
of very transparent contacts, it has been observed that 
carbon nanotubes may develop superconducting correla- 
tions, inherited from superconducting electrodes (prox- 
imity effect )& as well as intrinsic to large assemblies of 
massive ropesii*^. 

The superconducting instability is anyway in compe- 
tition with the so-called Peierls (or charge-density-wave) 
instability, which may induce a metallic-semiconducting 
transition caused by a lattice distorskm. The mean- 
field temperature Tp of such a transition has been es- 
timated by means of detailed calculations, and it is pre- 
dicted to increase as the radius of the nanotubes becomes 
smaller—. For tubes of typical radius, calculations find 
a very low (undetectable) value of Tp, while for thin- 
ner nanotubes it is predicted to be significantly larger 
and competing with the superconducting critical tem- 



In this paper we investigate the low-energy proper- 
ties of the (3,3) nanotubes by focusing on the instabil- 
ities of the Luttinger liquid behavior. We study care- 
fully the competition between the effective e-e interac- 
tion mediated by phonons and the Coulomb repulsion. 
The bosonization technique is applied in order to deal 
exactly with the interactions in the forward-scattering 
channels, while renormalization group methods are used 
to approach the low-energy instabilities of the system, 
driven by the backscattering and Umklapp interactions. 
Moreover, we pay also special attention to the experi- 
mental conditions reported in Ref. Il4l which lead to 
large arrays of nanotubes embedded in a zeolite matrix. 
This gives rise to a large screening of the Coulomb po- 
tential, which has no counterpart in the case of single 
nanotubesiSiiS. We study this effect by means of a gen- 
eralized RPA approach, showing that the long-range in- 
tertube coupling produces an efficient screening of the 
intratube interactions with small momentum-transfer. 

The most important result of the present study is the 
finding of two different low-energy phases characterized 
by critical (nonanalytic) behavior of the physical observ- 
ables. Under the conditions corresponding to the exper- 
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imental samples described in Ref. fill the critical be- 
havior is related to the vanishing of one of the Luttinger 
liquid parameters, and it is qualitatively consistent with 
the large diamagnetic signal observed in Ref. [lj. We 
also observe that this kind of singularity occurs well be- 
fore the development of any sizeable superconducting or 
charge-density-wave correlations in the electron system. 

We start by paying attention to the interactions medi- 
ated by the Coulomb potential, which provides a strong 
source of repulsion between electrons in single nanotubes. 
In the tubular nanotube geometry, the Coulomb poten- 
tial is given by£ 



yj {x - x') 2 + 4R 2 sin 2 [{y - y') 2 /2R 2 ] + a 2 

(1) 

where a z ~ 1.6 A and R is the nanotube radius. The 
screening by the environment of external charges is in 
general encoded in the dielectric constant k. For the 
sake of studying the nanotube transport properties, it is 
more convenient do deal with the one-dimensional (ID) 
projection of the potential onto the longitudinal dimen- 
sion of the nanotube. This is achieved by integration of 
the circular coordinate, upon which we obtain the ID 
Coulomb potential Vc(k) depending on the longitudinal 
momentum-transfer 

V c{k ) «^l g(^) , (2) 

k c is in general of the order of the inverse of the nanotube 
radius R, as it is the memory that the ID projection keeps 
of the finite transverse size. 

The relative strength of the Coulomb interaction is 
given by the dimcnsionless ratio between e 2 and the 
Fermi velocity vp, which for the (3,3) nanotubes is ap- 
proximately e 2 /vp ~ 2.9. This means that the Coloumb 
potential should give the dominant interaction in the 
forward-scattering channels, at least in single nanotubes. 
The processes can be classified depending on the low- 
energy linear branches involved in the electron-electron 
scattering 2 !. We recall at this point that the low-energy 
modes of the (3,3) nanotubes lie in a bonding and an an- 
tibonding subband that cross at two Fermi points (in the 
undoped system) with opposite longitudinal momenta 2 ^. 
We can distinguish in particular four forward-scattering 
channels, labelled by their respective couplings as repre- 
sented in Figs. ^and[3 In these processes there is a nom- 
inally strong Coulomb repulsion between the electrons, 
as they scatter without change of their chirality and the 
interaction strength is simply given by the Coulomb po- 
tential ©■ 

The interaction channels represented in Figs. n an( i|21 
are only part of the complete catalogue of scattering pro- 
cesses, that may be classified according to the chiralities 
and momentum-transfer for the interacting electrons' 2 -. 
It has become usual to assign respective coupling con- 
(7) 

stants gf to the interaction channels, in such a way that 




(b) 

FIG. 1: Small momentum-transfer processes corresponding to 
the couplings and g£ 




(b) 

FIG. 2: Small momentum-transfer processes corresponding to 
the couplings g^P and g^ . 

the lower index discerns whether the interacting particles 
shift from one Fermi point to the other (i = 1), remain 
at different Fermi points (i — 2), or they interact near 
the same Fermi point (i = 4). The upper label follows 
the same rule to classify the different combinations of 
left-movers and right-movers, including the possibility of 
having Umklapp processes (j = 3). As we will see, the 
couplings for the channels with large momentum-transfer 
2k f or change of chirality of the interacting electrons 
have an strength which is sensibly smaller than that of 
the forward-scattering couplings represented in Figs. ^ 
and El 
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FIG. 3: Schematic representation of the low-energy branches 
arising from the bonding and the antibonding subband of the 
armchair nanotubes. 



A very important point is that the system with 
just the forward-scattering interactions g^t\g^\g^ 

(2) 

and g\ ' is exactly solvable by means of bosonization 
techniquesSi22i2£. Thus, no matter that the Coulomb 
repulsion may place these interaction channels in the 
strong-coupling regime, the availability of an exact so- 
lution makes possible to capture the nonperturbative 
effects coming from the Coulomb interaction. The 
bosonization techniques make use of the fact that the 
forward-scattering interactions can be written in terms of 
the electron density operators corresponding to the differ- 
ent electron fields ^> r ia for the linear branches shown in 
Fig. |3 We adopt a notation in which the index r = L. R 
is used to label the left- or right-moving character of the 
linear branch, and the index i = 1 , 2 to label the Fermi 
point. The index a stands for the two different spin pro- 
jections. We may actually introduce the charge and spin 
density operators 

Pri{x) = -^=(^l ir (x)^(x) + ^ ril (x)^ ril ( X )) (3) 
Oriix) = -^=(^{x)^ rll {x)~^l tl {x)^ rll { X )) (4) 

As long as the Coulomb interaction and the interaction 
mediated by exchange of phonons do not depend on the 
spin of the interacting electrons, we will carry out the 
subsequent discussion in terms of the charge density op- 
erators p r i(x). 

It is convenient, for instance, to introduce operators 
for the sum and the difference of charge densities in the 
bonding and the antibonding subbands of the armchair 
nanotubes: 

PR±(k) = -j=(pm(k) ± p R2 (k)) (5) 
PL±(k) = ^=(p L2 (k) ± PL1 (k)) . (6) 

In terms of these operators, the hamiltonian for the 
forward-scattering interactions can be written in the form 

2 r^c 

Hfs = r«F / dk 2J ^ Pri(k)pri(-k) 

1 J - k ° r=L,Ri=± 



+p L+ (k) ( 3 f + <?f) p L+ (-k) 



+PR-(k) (fff - 9?>) PR-(-k) 
+ PL -(k) (gt ] - g ( 2 i} ) PL-(-k) 
+2 PR+ (k) (g { 2 2) + g[ 2) ) PL+(-k) 



-2pi 



-(*0 (9 ( 2 ] 



(2) 
94 



)PL-{-k)) 



(7) 



where k c stands again for the momentum cutoff dictated 
by the transverse size of the nanotube. 

We observe that the symmetric and the antisymmet- 
ric combination of the charge operators in the two low- 
energy subbands decouple in the hamiltonian (JJJ. This 
can be completely diagonalized by first introducing boson 
fields and $_(x) (and their respective conjugate 

momenta, 11+ (x) and EL (a;)) defined by 

d x <$>+{x) = ^(p L+ (x) + p R+ (x)) (8) 

«9 x $_(x) = V^(pl-(x) + p R -{x)) (9) 

The hamiltoninan can be rewritten then in the form 



Hfs — 



1 



dx 



K + (IL + (x)) 2 + 1 ^(d x $ + (x)) 2 



dx K_(n^(x)) 2 



^( dx ^(x)^) 



with renormalized velocities tt+, u_ and charge stiffnesses 
K +1 given by the expressions 

u±K± 



v F 



(sfi.gf)) (ID 
+ {g ( 2 2) ±gf ] )) (12) 

Upon application of the canonical transformation 



u ± /K ± = v F + ±(g^±g^ 




n 4 



=n. 



=n_ 



(13) 



the hamiltonian (|1C)|> becomes that of a couple of non- 
interacting boson fields. This picture, in which the low- 
energy excitations are just given by charge fluctuations 
(and spin fluctuations with unrenormalized velocity), is 
what characterizes the Luttinger liquid regime of the elec- 
tron systemSiSi. The charge stiffnesses and the renor- 
malized velocities in the two different charge sectors dic- 
tate for instance the thermodynamic and transport prop- 
erties, encoded into the compressibilities k±, the Drude 
weights D± and the dependence of the specific heat C v 
on the temperature 

« = ^ (14) 

7T U± 

2u±K± (15) 
1 1 



T 



3 u+ 



+ —) 



(16) 



The above picture will be complemented later by in- 
cluding the effect of the backscattering and Umklapp in- 
teractions, that tend to destabilize the Luttinger liquid 
regime at low temperatures. At this point, we stress that 
the robustness of such a regime depends on the strength 
of the Coulomb interaction. The consideration of the 
screening effects induced by the environment becomes 
then quite relevant, specially in the case of a large assem- 
bly of nanotubesiSiiS. For the nanotube array described 
m Ref.|Il the long-range Coulomb repulsion gives rise to 
a nonncgligible interaction between electronic currents in 
different nanotubes. If we label these by their respective 
positions 1 and 1' in the transverse section of the array, 
the intertube Coulomb potential can be expressed as 



Vl ,,,(*)«— if (|l-l'|fc), 

K 



(17) 



where k denotes the longitudinal momentum-transfer. 
We recall that the Bessel function K (x) is logarithmi- 
cally divergent in the limit x — ► 0. For 1 = 1', there 
is implicit a short-distance cutoff given by the radius 
of the nanotube, which leads then to the potential 
The important point is that the intertube Coulomb po- 
tential gives rise to quite significant screening effects at 
small momentum-transfer, which modify appreciably the 
strength of the forward-scattering couplings described 
above. 

In order to take into account the interaction among all 
the nanotubes in the array, we can adopt a generaliza- 
tion of the RPA scheme used in Ref. |27] for the study of 
2D layered sytems. In our case, the screened Coulomb 
potential V[J/(fe) has to satisfy the self-consistent dia- 
grammatic equation shown in Fig. 01 with 1" running 
over all the positions of the nanotubes in the array. We 
have then 



Vi,i'(fe) + n^vi,i. 



(18) 



where II stands for the polarization of each ID elec- 
tron system. This function is known exactly at small 
momentum-transfer* 2 *, and here it appears multiplied by 
the number of subbands n — 2 contributing at low- 
energies: 



IL(k,w) = 2n 



1 v F k 2 
ir uj 2 — dIF 



(19) 



Equation can be easily solved by introducing the 
Fourier transform of the Coulomb potential with respect 
to the nanotube position 1 in the 2D transverse section 
of the array. We define for instance 4>(p, k) by 



k) e l 



p.(i-i') 



(20) 



where d stands for the nanotube separation and BZ de- 
notes that the integration is over the Brillouin zone for 



i, V (r) ,v is V //' 




i\ v „, V 
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FIG. 4: Self-consistent diagrammatic equation representing 
the screening of the Coulomb potential between charges at 
nanotubes 1 and 1'. The sum over 1" accounts for the polar- 
ization effects from all the nanotubes in the array. 




FIG. 5: Dependence of the screened intratube Coulomb po- 
tential (in units of vf) on the dielectric constant k. 



the nanotube lattice in the transverse section of the array. 
We define similarly (j>( r \p, k) as the Fourier transform of 

Vil) (k). The expression of equation l|18(l becomes then 



m momentum space 

* (r) (p,fc) = tf(p,fc) 



Uc/>(p,k) 4> {r) (p,k) 



(21) 



In what follows, we will stick to the solution of this equa- 
tion at u> = 0, for the sake of giving a simpler description 
of the screening effects in the static limit. 

Within our RPA scheme, the screened Coulomb poten- 
tial becomes finally 



d~p 



(p,*) 



BZ 



1 + 2n0(p, k)/nVF 



ip-(i-i') 



(22) 

The most important property of V^) {k) is that it sat- 
urates at a finite value in the limit k — > 0. This is a 
reflection of the fact that, for distances much larger than 
the nanotube separation, the array screens effectively the 
Coulomb interaction as a 3D system. Thus, taking a 
value of the bare coupling e 2 w 2.9vp, we find that the 
intratube potential at vanishing momentum-transfer is 

V^\k = 0) w 0.75w F (for K = 1). This value depends 
slightly on the dielectric constant, as represented in Fig. 
We will take this strength of the screened intratube 
Coulomb potential at k = as an input for the values 
of the forward-scattering couplings , , and g^ 
within each nanotube in the array. 

A closer look at the above analysis shows that the 
Coulomb potential between nearest-neighbor nanotubes 
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FIG. 6: Dependence of the screened intertube Coulomb po- 
tential (in units of vf) on the dielectric constant k. 
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also gives rise to relevant interaction channels. The new 
couplings needed for a consistent description of the inter- 
action processes are catalogued in Fig. [7| They adhere 
to the same rules used to label the intratube couplings, 
but with a tilde to distinguish their intertube character. 
In particular, the intertube forward-scattering interac- 

(2) (2) 

tions g 2 an d 54 (preserving the chirality of the elec- 
trons within each nanotube) are affected by the screening 
effects described above, as they are corrected by the po- 
larization of the different nanotubes in the array. Their 
strength is given by the intertube Coulomb potential, 
which keeps only a nonnegligible value between nearest- 
neighbor tubes, after screening by the nanotube array. 
Taking again e 2 ~ 2.9vp, the intertube potential for 
|1 — l'| = d has a value at vanishing momentum-transfer 

V$(k = 0) « 0.007^ (for k = 1). The dependence 
of this strength on the dielectric constant is shown in 
Fig. El Although the relative strength of the intertube 
forward-scattering couplings g^ 1 and may appear 
small, these are however significant as they influence the 
scaling of the rest of interactions at low energies, as we 
will see in what follows. 

On the other hand, the intertube chirality-breaking in- 
teractions g^\g^\g^ and g^ give rise to a different 
kind of screening processes, of the type represented in 
Fig. [S] The effect of these processes cannot be cap- 
tured in the RPA scheme described above, as the po- 
larization with a change of chirality in the particle-hole 
pair diverges logarithmically at low energies. As long 
as the corrections depend on the ID cutoff, they give 
rise instead to new screening contributions to the scaling 
equations for the rest of interactions. We recall that the 
scaling equations for the interactions of single armchair 
nanotubes have been obtained in Ref. |^ (up to terms 
quadratic in the couplings). We have found that, after 
incorporating the corrections from the intertube interac- 
tions, the scaling equations for the intratube couplings 
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FIG. 7: Intertube interactions arising from the coupling be- 
tween electronic currents in nearest-neighbor nanotubes 1 and 
1' of a 3D array. The full (dashed) lines represent the prop- 
agation of electrons with right (left) chirality, and the labels 
1, 2 denote the respective Fermi points. The dashed lines 
(without arrow) stand for the intertube Coulomb potential. 
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FIG. 8: Second-order process renormalizing backscattering 
interactions at a given nanotube 1 through the coupling with 
the nearest-neighbors 1' in a 3D array of nanotubes. The 
dashed lines (without arrow) stand for the intertube Coulomb 
potential. 
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The variable Z stands for minus the logarithm of the en- 
ergy (temperature) scale measured in units of the high- 
energy cutoff of the ID model E c ~ vpk c (of the order of 
~ 0.1 eV). The large coefficients in front of the intertube 
contributions arise from the number of nearest-neighbors 
of each nanotube in the 3D array. Furthermore, we have 
also incorporated a nonperturbative improvement of the 
equations by writing the exact dependence of the scal- 
ing dimensions on the forward-scattering couplings, ex- 
pressed in terms of the K + , K_ parameters. 

The new intertube interactions are themselves cor- 
rected by processes that diverge logarithmically at low 
energy, and that give rise to respective contributions to 
the scaling of the intertube couplings. We can focus on 

the analysis of (jf^si^isi an d si 3 \ which have values 
given by the intertube potential (|17|l at the initial stage 
of the low-energy scaling. As long as this potential decays 
exponentially for |1 — l'\k > 1, we can neglect the influ- 



ence of other intertube interactions with 2\zf momentum- 
transfer. The second-order diagrams that renormalizc 
the above intertube couplings consist of particle-particle 
processes or particle-hole loops involving a change of chi- 
rality, as illustrated in Figs. l9landlTUl Some of the contri- 
butions depend on and g\ . This means that these 
couplings have to be taken into account for a consistent 
description of the low-energy scaling. The complete set 
of scaling equations for the intertube couplings becomes: 
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Following the flow of the scaling equations, the 
backscattering and Umklapp interactions are progres- 
sively enhanced as the theory is scaled down to low en- 
ergies. At the initial stage of the renormalization, the 
values of the couplings are dictated by the Coulomb in- 
teraction and, in the case of intratube couplings, also by 
the effective interaction mediated by phonon-exchange. 
Regarding the Coulomb interaction, its contribution to 
forward-scattering couplings is given by the RPA cal- 



culation exposed above. Thus we have g\ 
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0)/vf, following the trend 

shown in Fig. while gf ] = ~gf ] = V$(k = 0)/v F , 
with the potential for nearest-neighbor 1, 1' represented 
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FIG. 9: Second-order diagrams with logarithmic dependence 
on the frequency renormalizing the intertube interac- 
tion. The wavy lines stand for intratube interactions and the 
dashed lines (without arrow) for interactions between nearest- 
neighbor nanotubes 1, 1' in a 3D array. The rest of the lines 
keep the same meaning as in Fig. |7] 



FIG. 10: Second-order diagrams with logarithmic dependence 

(3) 

on the frequency renormalizing the intertube §2 interac- 
tion. The wavy lines stand for intratube interactions and the 
dashed lines (without arrow) for interactions between nearest- 
neighbor nanotubes 1, 1' in a 3D array. 



in Fig. For the chirality-breaking processes, the dif- 
ferent symmetry of ingoing and outgoing electron modes 
implies also a significant reduction of the Coulomb po- 
tential, as evaluated in Ref. 0- The result is that, for 
the armchair (3,3) nanotubes, there is a repulsive com- 
ponent in the backscattering and Umklapp interactions 
at small momentum-transfer whose strength can be es- 
timated as w 0.23«f/k. This applies to the couplings 

<?2 > 9^ > 9^P an d ffi j as wen as to tne corresponding in- 
tertube couplings at small momentum-transfer. For the 
intratube couplings with the large momentum-transfer 
2kp, the contribution by the Coulomb interaction can 



be obtained from the Fourier transform of the potential 
JTJ. The strength estimated in this way turns out to be 
« 0.06w f /k. 

Dealing now with the contribution from the effective 
phonon-mediated interaction, we note that the situation 
is reversed, and that the backscattering and Umklapp 
couplings at vanishing momentum-transfer get a smaller 
component than the couplings for momentum-transfer 
around 2kp. Previous estimates had already found that 
the ratio between these two different strengths is approx- 
imately l/3*&. More recent calculations of the phonon 
spectrum by means of density functional theory have 
led to a similar proportion 1 -^. More precisely, it has 
been found that the contributions of all the phonons 
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K 

FIG. 11: Phase diagram showing the different low-energy 
instabilities in the array of (3,3) nanotubes, depending on 
the dielectric constant k and the effective coupling A of the 
phonon-mediated interaction. The three thin curves corre- 
spond to constant l c lines in the phase diagram: l c = 4, 6, 8 
from top to bottom. 



with momentum-transfer 2kp add to an effective cou- 
pling A w 0.1, while the contributions of the phonons near 
the zone center give an effective coupling w A/ 3. These 
more accurate estimates are about three times smaller 
than those quoted in Ref. 28. Anyhow, we will con- 
sider the phase diagram of the (3,3) nanotubes by span- 
ning a suitable range in the scale A of the two effective 
couplings, covering the values between the different esti- 
mates in Refs. ^| andl28l 

In order to determine the electronic instabilities that 
may appear at low energies, we have solved the set of 
scaling equations P% |l -l|3T jl . taking initial values for the 
couplings according to the above discussion. The cou- 
plings approach in general a regime where they grow 
large as I — > oo. Regarding the forward-scattering in- 
teractions, g\ becomes increasingly repulsive, leading 
to a singularity characterized by either the vanishing of 
K-y or the divergence of K_ , depending on the region of 
the phase diagram. The scaling of the interactions stops 
at the low-energy scale uj c = E c e~ lc corresponding to the 
point l c where the singularity is reached. This has ac- 
tually the character of a critical point, since it gives rise 
to the opening of a branch cut and nonanalytic behavior 
in the corresponding Luttinger liquid parameter K + or 
K_—. We have plotted in Fig. ^]the phase diagram of 
the (3,3) nanotubes showing the two different regions of 
singular behavior, as a function of the dielectric constant 
k and the effective coupling A of the phonon-mediated 
e-e interaction. 

In general, the enhancement of backscattering and 
Umklapp interactions upon scaling may lead to a large 
growth of electron correlations, pointing at the tendency 
towards long-range order in the electron system. We have 
analyzed this possibility through the computation of dif- 



ferent response functions 

/-f- OO f-h 
dt / dxe iut e ikx (TO(x,t)O(Q,0y) 
-oo JO 

where the pair field O characterizes a particular type of 
ordering. In the system under consideration, the most 
important correlation functions are found to be given by 
the following (Fourier transformed) fields: 

DWtlt (k « 2k F ) = -±= [*Ju«(P - k)a a ^ L2fj (p) 

p,<x,(3 

DW ,Ak « 0) = [*ka(P- k)a^ Lip (p) 

+ *k«(p-*K ,/J *M/j(p)" , 

OscAk « 0) = a [*«ia(-P + k)a- a -^ L2p (p) 

P,ct,l3 

+ ^R2a(-p+k)a- a ^ Lll3 {p)] , 

where, for density wave (DW) operators, fi = corre- 
sponds to a charge-density wave (CDW) and fi = 1,2,3 
to a spin-density wave (SDW); while, for superconducting 
(SC) operators, fi = stands for singlet superconductiv- 
ity and fi = 1, 2, 3 for triplet superconductivity (cr^'' 3 are 

the Pauli matrices, with a^'^ = 12x2)- 

The derivatives with respect to the frequency of the re- 
sponse functions obey actually scaling equations^, whose 
solution allows to compare the relative growth of the dif- 
ferent electron correlations. By looking at the low-energy 
scaling, we have checked however that none of the re- 
sponse functions shows a very large growth, down to the 
point where the scaling flow breaks down due to the sin- 
gularity in the Luttinger liquid parameter. This singular 
behavior occurs therefore before the appearance of any 
tendency to long-range order in the electron system. This 
is illustrated for a typical instance in Fig. ^1 where it 
can be observed that only the CDW correlations with 
vanishing momentum show a significant growth as the 
critical point is approached. 

The (3,3) carbon nanotubes may fall therefore into two 
different low-temperature phases, whose physical proper- 
ties are dictated by the vanishing of K + and the diver- 
gence of K_ respectively We remind in particular that, 
for the experimental setup of Ref. 113. a reasonable choice 
of the parameters is k ~ 2 -j- 4 and A ~ 0.1, correspond- 
ing to the K + = phase. In this case the temperature of 
transition to the new phase results strongly dependent on 
the dielectric constant of the environment, ranging from 
T c - 10" 1 K (at k w 2) to T c - 10 K (at k ps 4). 

The behavior of the response functions follows in gen- 
eral the trend shown in Fig. 1121 and, while the density- 
wave correlations tend to grow by approaching the crit- 
ical value l c , the superconducting correlations remain 
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FIG. 12: Plot of different response functions at k = 2 and 
A = 0.1. 

small anyhow. This finding seems to rule out the possi- 
bility of having superconducting correlations in the (3,3) 
nanotubes, at least under the physical conditions con- 
sidered in the present paper. We coincide in this re- 
spect with the conclusions reached in previous analyses 
by means of other computational methods 15 - ^ We have 
found however that the appearance of a CDW instabil- 
ity is precluded by the breakdown of the Luttinger liquid 
behavior, which cuts off the growth of the different cor- 
relations. The destabilization of the Luttinger liquid is 
favored in this respect by the screening effects from the 
3D array of nanotubes, which are responsible of bringing 
the Coulomb interaction into the weak-coupling regime. 

Our results confront the claim that the experimen- 
tal signatures reported in Ref. Q should provide ev- 
idence for a superconducting transition in the small- 
diameter nanotubes. This interpretation has been also 
challenged by studies of the electron correlations in the 
(5,0) nanotubesi^iSi. We have shown that, even consid- 
ering the large screening effects from the arrays of nan- 
otubes in the experimental samples, the effective e-e at- 
traction arising from the exchange of phonons is not large 
enough to support the appearance of superconducting 
correlations in the (3,3) nanotubes. If the Coulomb inter- 
action is further screened by a suitable variation of the di- 
electric constant of the medium, the system is driven then 



into the phase characterized by the divergence of K_ and 
the related compressibility k_ , as shown in Fig. ^] This 
has the same character that the instability given by the 
Wentzcl-Bardccn singularity, where the divergent com- 
pressibility is the signal of the spatial separation of the 
system into regions with different electron density 31,32 . 

Anyhow, the experimental conditions of the samples 
described in Ref. [Ij seem to place the system in the re- 
gion of the phase diagram characterized by the vanishing 
of K + . This leads to the vanishing of the conductivity 
at the point of the transition, as follows from Eq. JTSJ. 
Moreover, it also gives rise to the vanishing of the tun- 
neling conductance, which is directly related to the tun- 
neling density of states n(e). Within the Luttinger liquid 
framework, the latter follows the low-energy behavior 

n(e) ~ e (K++i/K + +K-+i/K— 4)/8 (39) 

The depletion of the density of states given by Eq. i|39|) 
at vanishing K+ is consistent with the appearance of the 
pseudogap observed experimentally in the measures of 
the I-V characteristics reported in Ref. The critical 
point characterized by the vanishing of K + does not de- 
scribe however a conventional metal-insulator transition, 
as long as the compressibility given by Eq. I|14l) remains 
finite at the point of the transition. As analyzed in Ref. 
13.1 the critical point implies actually a divergent diamag- 
netic susceptibility, as a consequence of the development 
of very soft modes in the sector of electron current exci- 
tations. Therefore, the phenomenology derived from the 
K + = phase of the (3,3) carbon nanotubes seems to 
be consistent, at least qualitatively, with the main ex- 
perimental signatures reported in Ref. Il4i Further ex- 
perimental work would be needed to clarify the existence 
of such a critical point in the (3,3) nanotubes, its physi- 
cal properties, and its stability under changes of relevant 
experimental conditions. 
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